Ultralocal Energy Density in Massive Gravity 

Vladimir O. Soloviev and Margarita V. Tchichikina 
June 29, 2011 

Abstract 

We provide a space-time covariant Hamiltonian treatment for a finite- 
range gravitational theory. The Kuchai^ approach is used to demonstrate 
the bimetric picture of space-time in its most transparent form. This 
Hamiltonian formalism is applied for the straightforward realization of the 
Poincare algebra in Dirac brackets. It uncovers the simplest form of the 
Poincare generators expressed as spatial integrals of ultralocal quantities 
constructed pure algebraically by means of the two space-time metrics. 

1 Introduction 

For many reasons, both observational and motivated by pure theory, modifica- 
tions of the gravitational theory attract a lot of attention now. In particular, we 
mean those changes in the standard model of gravity that modify its infrared 
behavior. It is possible to achieve the transformation of the massless theory into 
the massive one in different ways. The discussion may be found, for example, 
in articles [U [2 [3l |4j [5l |6] . In this paper we deal with a variant of gravitational 
theory which is historically connected with works [71 HI |9] . 

On the ground of the fact that the linearized kinetic energy contains a neg- 
ative mode [TU] sometimes it is said that this model does not deserve futher 
study. But as some other authors [TTJ [TU [T31 [H] we can not agree with this 
statement. It is necessary to study the behaviour of solutions not only in linear 
case, but in the full theory, as it can be rather different. 

The content of the theory becomes in many aspects more transparent after 
expression in the general language of Hamiltonian formalism. With this pur- 
pose we demonstrate here the construction and elementary applications of the 
canonical formulation for the Relativistic Theory of Gravitation (RTG) |TS]. We 
stress that our main interest is in the generally covariant canonical formalism 
developed by Kuchai^ [W\ and not in the ADM variant [T7]. An essential dif- 
ference of Kuchai^'s approach is that four components of (pseudo)Riemannian 
metric tensor and four lapse and shift functions arc not identified and play dif- 
ferent roles. Precisely, the momenta conjugate to the four metric components 
are zero and that leads to four primary constraints, whereas those metric com- 
ponents are persent in the secondary constraints. All the four pairs of variables 
are excluded after solving the second class constraints, whereas lapse and shift 
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functions stay in the Hamiltonian, providing its general coordinate invariance 
and freedom to choose the foHation of spacetime by spatial hypersurfaces. Surely 
both the Kuchaf approach and ADM one are based on the fundamental Dirac 
results, presented in brilliant form in his Yeshiva Lectures (TB]- Here we apply, 
in fact, just in the same order as in the book [TS], both the variant called the 
quantization on curved surfaces and the other one called the quantization on 
flat surfaces. 

The bimetric view of General Relativity (GR) first appeared in publica- 
tion [7]. After expressing a bimetric theory in the Hamiltonian form we im- 
mediately encounter with the causality issue. In the gravitational theory |15) 
discussed here this problem is solved by introducing an extra causality pos- 
tulate. This postulate requires that the null cone defined by the dynamical 
(pseudo)Riemannian metric should be inside or coincide with the null cone de- 
fined by the flat background metric. Then any hypersurface which is spacelike 
according to the flat metric is also spacelike in the (pseudo)Riemannian metric. 
As covariant Hamiltonian approach allows a free choice of time one always can 
take it so that its direction will be timelike in both metrics. We demonstrate 
that the total energy density is proportional to the square of the unit normal 
to the constant time hypersurface. Of course, the unit normal is defined on the 
base of one metric and it is squared with another one. So the time problem and 
energy problem become united. 

At last, most of physicists believe that the gravitational field as other funda- 
mental fields subjects to quantization. The difhculties of metric tensor quanti- 
zation in General Relativity are well known. Perspectives of quantization of the 
alternative gravitation theories deserve further study. Canonical quantization is 
the first in the list of methods developed to unite quantum and classical theory, 
both historically and practically, so we start here with a discussion of the Hamil- 
tonian formalism also having in mind a program of quantization. In Ref. [19] the 
bimetric gravity has been treated as a zero graviton mass theory, but later that 
version has been declined. The new statement of this problem, as will be seen 
below, leads to rather different results. There are no first class constraints in the 
theory now, the number of degrees of freedom increases and Poincare invariance 
results in ten non-trivial integrals of motion. Some preliminary presentation of 
this work appeared in Refs. [20l |2T] 

In Section 2 we introduce our notations and the Lagrangian for the con- 
sidered theory of gravitation. Then we apply to all tensors and vectors the 
spacetime covariant 3 -I- 1-decomposition developed by Kuchai^ [16j . Spacetime 
is treated as foliated by one-parametric family of spacelike hypersurfaces. The 
time here is not a coordinate, but a parameter continuously numbering hyper- 
surfaces. This formalism exploits simultaneously two coordinate systems: the 
first is general and arbitrary, the second is induced by the foliation. General 
spacetime coordinate invariance so is respected. The pecularity of bimetric 
picture is a strong reqirement that hypersurfaces are to be spacelike in both 
metrics. There are also two bases for decomposition of tensors corresponding 
to the two spacetime metrics. These decompositions are applied both to proper 
spacetime tensors and to tensors depending on the foliation. Practically, more 
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attention is paid to the decomposition based on fiat background metric, though 
sometimes it is rather suitable to use another basis. 

In Section 3 (3 + l)-decomposed tensors are substituted into the Lagrangian 
and after this transformation of variables the Hamiltonian formalism is con- 
structed. The Lagrangian occurs to be a sum of two terms, one of them coin- 
sides with the term of General Relativity, the other is constructed from the two 
metric by purely algebraic means. So, the transformation of variables is not 
difficult. In contrast to one-metric theories, such as General Relativity, all ten 
components of (pseudo)Riemannian metric enter the Lagrangian and should be 
treated as dynamical variables. But velocities of four of them do not enter the 
Lagrangian and so we get four primary constraints. Conservation of them leads 
to four secondary constraints and all eight constraints occur second class. It is 
easy to solve four secondary constraints and express four variables as functions 
of other variables. After construction of Dirac brackets all four pairs of above 
mentioned variables and all constraints may be excluded. In the same time 
Dirac brackets for the rest gravitational variables coincide with their Poisson 
brackets in General Relativity. The Hamiltonian of the gravitational field has 
a standard for field theories on fixed background form, i.e. it is linear in lapse 
and shift hmctions. 

In Section 4 we take the scalar field as a representative of matter. The 3-1-1- 
decomposition again leads us to the field Hamiltonian linear in lapse and shift 
functions. The gravitational interaction(minimal interaction) of scalar field is 
provided by (pseudo)Riemannian metric. Only uniting both Lagrangians we 
obtain the consistent interaction model. We present the full system of canonical 
equations generated by the Hamiltonian and Dirac brackets. 

Section 5 is devoted to the Poincare group and to generation of its transfor- 
mations by the Hamiltonian. As Hamiltonian depends on four functions N, N'^ 
which are not dynamical variables the evolution depends on their choice. Start- 
ing the evolution from a given hypersurface we are free to choose the following 
hypersurfaces (of course by preserving smoothness). The evolution parameter 
is not always possible to be treated as some kind of "time", in general, it is 
simply a parameter of some smooth coordinate transformation. In particular, 
the Hamiltonian can generate purely spatial coordinate transformations in a 
hypersurface if we put lapse function to zero. We start with a calculation of 
Dirac braket for two Hamiltonians with different lapse and shift functions. This 
bracket almost gives the celebrated "hypersurface deformation algebra" derived 
by Dirac. For the full correspondence the Hamiltonians should not contain other 
functions which are not canonical variables then lapse and shift. In our case such 
a function is a spatial metric induced on hypersurfaces by the flat background 
metric of spacetime. Theories with the background metric may be extended 
by introduction of four pairs of canonical variables: embedding variables and 
their conjugate momenta, then the exact Dirac formulas are valid. Also these 
formulas are valid in General Relativity because it does not contain any back- 
ground metric. In our case the algebra appears when we restrict ourselves by 
transformations which are motions for the background metric (Killing vectors) . 
For the flat background they are transformations from the Poincare group. In 
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order to preserve not only spacetime background metric but also spatial metric 
induced by it on hypersurfaces it is necessary to restrict formalism by those 
which are flat in background metric. And the last restriction is not necessary - 
for simplicity we allow only Cartesian coordinates on the hypersurfaces. 

If we put a requirement that the empty flat space has exactly zero energy 
density, then the Poincare algebra will have a central extension and central 
charges. In other case we will have an algebra without central charges but 
the empty Minkowski spacetime will have some small constant positive energy 
density. 

We should like to draw attention to the fact that energy and momentum 
densities are ultralocal here and this gives a straightforward connection between 
the sign of total energy density and the bimetric causal structure of spacetime: 
energy density occurs proportional to the square of unit normal to the hypersur- 
face. Of course, this unit normal should be defined by means of background flat 
metric and its square is calculated by means of another, (pseudo)Riemannian 
spacetime metric. The reasonable choice of foliation, therefore, should provide 
us with positivity of the total energy density for gravitational field and matter. 

2 The RTG Lagrangian, 3+1-decomposition of 
tensors and a new form of the Lagrangian den- 
sity 

We start with the RTG Lagrangian [TS], which allows in more or less stan- 
dard way to derive the Hamiltonian and the Poisson brackets. The RTG 
Lagrangian contains nondynamical flat metric /i^^ together with dynamical 
pseudo-Riemannian metric g^^, this is the main difference from the GRT. Tak- 
ing the velocity of light equal to unity, we can take the Lagrangian density of 
the gravitational field following Ref. [15] as 

where dots denote surface terms (4- divergences) , Greek indices take values from 
to 3, G is the gravitational constant, g = det((7^[^), g^^^ = ^J—gg^^, h = 
det(/i^^), R is the space-time scalar curvature given by metric 17^1,, m is the 
graviton mass. The signature of space-time is (—1, 1, 1, 1). 

We can also present Lagrangian density ^ in an equivalent (up to diver- 
gences) form 




where 



^r^i/ = - = -g^" {D^g^^ + D^g^^ - D^gf,^) , (3) 
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f^^ are the ChristofFel symbols of the pseudo-Riemannian metric, F^^ are 
ChristofFels of the flat metric, and is the covariant derivative compatible 
with the flat metric. 

As independent variables in the action constructed on the base of Lagrangian 
densities ([ll or (l2|), which should be varied, one can take, for example, 10 com- 
ponents of pseudo-Riemannian metric 5^1/. To simplify the formal calculations 
the following tensor objects will be useful 



r ^ ^Z"- (4) 

In constructing the Hamiltonian formalism it is necessary to chose the arrow 
of evolution, i.e. the physical time direction. This does not mean violation of 
the general covariance and fixing one of space-time coordinates X" as time. It 
was firstly shown in Kuchaf 's works |16) how to make the canonical formalism 
space-time covariant. We will follow this approach here, and we will also use 
some results derived in paper j22) . 

The fixed instant of physical time corresponds to some spacelike hypersurface 

X" = e"(a;'), (5) 

where are independent coordinates on the hypersurface, Latin indices take 
values from 1 to 3. Unlike the case treated in publications [TB], we deal here with 
two metrics for the spacetime, not one, and we demand that our hypersurfaces 
are to be spacelike in both metrics. It means imposing two conditions valid 
at any point of the hypersurface 

J^j{x'')dx'dx^ > 0, rirj{x^)dx''dx^ > 0, (6) 

where two different induced metrics are involved: 

Itj = g^.•^e'^^eJ, Vij = h^^e'^.e'^^. (7) 

It is evident that in general case metric 77.^ , in contrast to metric /i^i/, is not 
flat. 

We will assume that the spacetime may be filled with such slices of constant 
physical time, i.e. presented as one-parametrical family of spacelike hypersur- 
faces: 

X"' = e°'{x\t), (8) 
and we can introduce some vector field 

iV" = ^, (9) 

which is timelike in both metrics of spacetime 

gapN^N^ < 0, KpN'^N^ < 0. (10) 



^This requirement is justified by the postulate of causaUty in the RTG. 



5 



In order to provide 3 + 1-decompositions of spacetime tensors we are to have 
a basis related to a fixed time hypersurface. For example, we may use for this 
purpose four spacetime vectors (A^",e"). But really we need two bases: one 
connected to metric g^^, and another connected to metric h^^. This doubling 
occurs when we go to lower indices 

Np^hp^N'^, Np=gp^N'', (11) 

and analogously, Cai and e-d appear. It is more suitable from the technical view 
to introduce two other bases, where unit normals to the hypersurface are taken 
as timelike vectors (n",e") and (fi",e"). They are defined by the following 
conditions 

Kpn'^e'^^^=Q, ha^pn^n^ = -I, (12) 

5a0n"e(l = 0, go/3n"n^ - -1. (13) 

Now we are able to apply 3+1-decomposition to different tensors, for exam- 
ple. 



f"" - f^^n'^n'' + f^hi''e':',+ f^ey + f^e'^^e'^.j, (14) 



where 



N^-n^N'', N'^e'^m, N^-n^m, N' = el^N^', (15) 



5^^ = n^n.5^^ 5^^' - .g^^ = -n^ei5^^ g^^ ^ e^eig^'' , . . . (16) 
It is easy to find the linear relations between the two bases 

V-9 

and correspondingly the linear relations between the components as 



^=-J^\fl^^ N' = N'^-jTlN. (18) 



3 Construction of the canonical formalism for 
the gravitational field 

In order to transform the Lagrangian density to the desired form ([T]) it is neces- 
sary to express in a new way the two terms, where the first one does not contain 
neither the graviton mass, nor the flat metric. Therefore for the first term it is 
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enough to apply the standard transformations [TB] used in GRT. Then up to 
surface terms we obtain 



N 



16nG f^^^ 



(19) 



where R is the scalar curvature of the hypersurface, derived by means of met- 
ric "fij, Kij is the second fundamental form of the hypersurface in pseudo- 

" ,KkiY'"-f'^- Taking 



Riemannian geometry defined by metric g^i, 
into account formula 



1 



(20) 



and after substitution of decompositions ([T^ , the second term takes the follow- 
ing form 



NVV' 



-1 - 



1 7 /I 



1 



(21) 



It is seen that this expression does not contain any velocities and so it is irrele- 
vant for definition of momenta. 

The surface terms in Lagrangian density (total derivatives and spatial diver- 
gences) do not influence the symplectic structure and so the Poisson brackets. 
Boundary conditions at the spatial infinity in their turn are taken so, that the 
pseudo-Riemannian metric tends to flat Minkowski metric, and the hypersur- 
faces tend to hyperplanes. So we arrive at the following form of the gravitational 
field action 



dt I (fx I - 



tl i?3 



N 



(i? - + S-pK^ 



+ 



2/J 



(22) 



Here ^ij{x'',t), f^-^{x'',t), /^*(x'^, t) are the independent variables to be varied. 

Flat spacetime metric /i^^ (X") and functions e"{x^,t) which parametrize the 
spacelike hypersurfaces are treated as given and so not to be varied. Base vectors 
n"(x% t), ef (a;*, i) and vector 7V"(x% t) are in their turn expressed through those 
functions. Quantities Kij{x^^t) are given by the formulas well-known from the 
canonical formalism of General Relativity 



1 
2N 



N, 



3\i 



(23) 



where functions N, TV* in their turn are expressed through N, iV' and /^^, 
by equations ^8]) . Here the vertical bar denotes the covariant derivative 
determined in the Riemannian geometry of 3-dimensional space by metric "/ij . 
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In order to avoid confusion of mixing momenta with the famous mathematical 
constant tt let us hide this constant in a new constant n = IGttG. 

Starting from action we obtain for the conjugate momenta the following 
expressions: 



7r± 



dC 
dC 

w 

dC ^ 



= 0, 
= 0, 



dC dKi, 



(24) 
(25) 
(26) 



It is evident from the above that equations p4|) and p5|) are primary constraints 
in Dirac's terminology [18], and so these equations are to be added with arbitrary 
multipliers to the Hamiltonian. Therefore we obtain 



ii3 



(27) 



where velocities should be expressed through moments as follows 



(28) 



Now, after finishing this procedure, the Hamiltonian takes the following form 
(up to surface terms) 



H = j d?x {NH + N"H, + A^TTi + AV,) , 

fl3 



(29) 



where 



n = - 



1 



nj J_ nj 



2 + 2/^^ f^^v\2 



1 7 



n 



-27rV, 



The canonical Poisson brackets given as follows 

5F 6G 6F SG SF SG 



{F,G} = J 

fl3 



Sjij Jtt^J (57r_L Sf^' Sir, 



(F^G) 



(30) 
(31) 
(32) 

(33) 
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let us to present the hamiltonian equations in a standard form: 



(34) 

(35) 
(36) 



Next we are to check whether the primary constraints (|24p , (1251) are compat- 
ible with the equations of motion, in order to guarantee that time derivatives 
7r±,o and tt^^o are also zero. 

As conjugate variables /-'--'- and enter the Hamiltonian algebraically we 
obtain secondary constraints in the form of algebraic equations: 



dn 



= 0, 



dn 



-0, 



gf±± ' Qf±t 
which can be solved in elementary way and they give relations: 



K 



(37) 



(38) 



K rj 



It is easily seen from the solved form of the secondary constraints that their Pois- 
son brackets with the primary constraints are nonzero, i.e. all the constraints 
are second class and so may be excluded by introducing the Dirac brackets. Here 
Dirac brackets may be obtained from Poisson brackets (|33p by simple exclusion 
of terms with variables {f-^^,n^_) and [f^'^^iTi) 



{F,G}d = Jd' 



6F SG SF 6G 



(40) 



After substituting the solutions of constraint equations into the hamiltonian we 
obtain 



H 



N 



In 



m^y/jy 7 / 1 



K rj 



(41) 



We arrived at the Hamiltonian depending on canonical variables 7i j , tt'^ , 
and on the prescribed metric rjij (determined from the fixed spacetime metric 
and the functions parametrizing hypersurfaces) . The constant provides vacuum 
energy normalizing, i.e. when the Riemannian metric coincides with the flat 
one g^^ = h^^, we have H = for any spacelike hypersurfaces. 
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Let us add that it is possible to do not exclude all constraints from the 
Hamiltonian and exploit the Dirac brackets from the very beginning. We also 
can at first make secondary constraints compatible with dynamics of equations 
- ([55]) - ((55)) , this will give us a chance to determine the Lagrangian multi- 
pliers through canonical variables, but this will not lead us to new constraints. 
After substituting these Lagrangian multipliers into equations (|55|) - we 
will obtain relations which are equivalent to the famous harmonicity conditions 

D^r-" = 0. (42) 

But the Dirac brackets exclude the independent role of these equations and 
these equations become a consequence of the Hamiltonian equations of motion 
generated by Hamiltonian (HT]) and brackets 



4 Scalar field as an example of gravity source 

The formalism developed above will be obviously incomplete without a demon- 
stration of extending it onto matter fields. Let us suppose that their interaction 
with gravity is minimal, then the matter fields Lagrangian density Cm will de- 
pend on the set of fields (f>^{X"), on their spacetime coordinate derivatives and 
on metric gf^^{X°')^ transforming as the scalar density under general coordinate 
transformations. The Hamiltonian form of matter fields action can be obtained 
through Kuchaf's procedure ^16) . and as a result we get 

*2 

Sm = J dtj d^x {TTAq^^o - N-Hm - N"HM^) ■ (43) 

The combination of this action with the gravitational action (1221) results simply 
in adding of quantities T-Lm and to previously obtained % and "Hi, which 
depend on gravitational variables only. 

To illustrate this we consider the scalar field with the Lagrangian density 



= -V-5 ( 7;9''''d^<pd,^ + U{<jy) ) (44) 



The transformation to 3-|-l-notations and the Legandre transform may be done 
separately of the gravitational contribution 

Cm - -N^^ (^/^^</>a0,± + 2/^^,^^. + (j^' f^' - ^) ^I^A., + 

(45) 

where 

0,0 (46) 
The momentum is determined as usual 

fir 

^'^^^ = VV - , (47) 
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and the velocity may be expressed through the momentum as fohows 



N ( ■ f-L' 



(48) 



After providing the corresponding Legandre transform the action for scalar field 
will take the form (j43]), where 



M 



V7 



(49) 



Therefore, after combining actions of the gravitational and the scalar fields we 
get the same primary constraints (^51 . as in the case of pure gravity, 

and the full Hamiltonian with constraints will have the same form (l29l) . The 
procedure of constraints exclusion also does not change and the final form of 
the Hamiltonian remains the same (|4ip. where now 



1 



1 



H = -= I --7i? + «(Sp7r^-— ) + ^ + -77'^a,(/.9,(/. + 7[/(0) 



i?3 



K - 2 ' 2 2 

SF SG SF SG SF SG SF 5G 



S^ij Sn'-^ S<j) Suff, S^ij Sn^ S(f> 



(50) 



It is easy to be convinced that the hamiltonian equations of motion for the 
system of interacting scalar and gravitational fields are of the following form (in 
the case U{<l)) = l/2AP(jP): 



5,0 



d'x(N{(j,,-H}D + N''{^,'Hk}D 



(51) 



(iVV77'''a,(^),, - Ny^M^(b + (^V^),„ 



d-^x 



(^N{j,,,n}D + N''Hj,nk}D 



2N . TT, 

'x(^{n'' ,m}D + N'^Itt^' ^nkjn) 



fl3 



i'' + i^m {-/''W - 7''7") 



(52) 



(53) 
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2 J.2 



1 - N ■ 



(54) 



It is evident that the difference with the corresponding equations of General 
Relativity appears in the last equation only and has the order of magnitude 
0{w? /k). The dependence of quantities /^^, on variables (|39t . whose 
Hamiltonian equations have been derived above, may be ignored in the calcu- 
lating of Dirac brackets acording to formulas p7p . 



5 Poincare group in Hamiltonian formalism of 
the RTG 

Among all variants of the hamiltonian evolution arising from arbitrariness of 
functions N{x), N'^{x) Hamiltonian (|4ip contains special transformations pre- 
serving the Minkowski metric. So, if we take hyperplanes as hypersurfaces and if 
we take Cartesian systems of coordinates on them, we will have on these hyper- 
surfaces metric r/ij induced by Minkowski metric of spacetime (7) in the simplest 
form r]ij = Sij , and we will have transformation functions (|15p as follows 

N = Akx'' + a, N' = A,kX^ + a\ (55) 

where 

Aik = ~Aki- 

Then Hamiltonian (|4T|). according to its linearity in functions N{x), N^{x), will 
take a form 

H ^ P°a~ P'a' + M'^Ak + ^AP'^A^k, (56) 



where 



,2 



m 

K 



f^'d^x = - J n,d^x, 
M''' = J {x'f^'' -x''f^')d^x = J {x'^m- x"Hk)d^x, 

M" = ~— I x\l + f^^)d'x. (57) 



K 



The meaning of these operators is clear as they are special cases of the general 
Hamiltonian and correspond to different choices of coordinate transformations 
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which they generate: P° refers to the time translation, - to the spatial trans- 
lations, M*'^ - to the spatial rotations and M'^ - to Lorentz boosts. Our nota- 
tions are chosen in order to make easy comparison with the analogous formulas 
from paper |23| , where Poincare algebra has been considered in the asymptoti- 
cally flat spacetime of General Relativity. 

But before reducing our Hamiltonian to such restricted form it is useful to 
derive the algebra of Dirac brackets (|50p for the general Hamiltonians. Let us 
take 



where we treat /^^, as functions having zero Dirac brackets. This is jus- 
tified by the well-known statement that second class constraints may be taken 
into account both before and after calculations of Dirac brackets. So, we will 
account for the first half of constraints ([M]). (I^5t before, and the second half 
([M)) - after. As usual all surface integrals in integration by parts are dis- 
carded. It is reasonable for island systems when radiation is allowed only inside 
of them, not at infinity where the pseudo-Riemannian metric tends to the fiat 
one with the rate given by the Yukawa behaviour. 

The results of calculations may be presented in the form suitable for com- 
parison with the analogous formulas of GR: 



or in the form corresponding to parametrized field theories at the background 
of flat metric: 



H 





(59) 




A' 



A = 



(60) 
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where Csfjij is the Lee derivative of metric rjij in the direction of vector field a. 
The difference from GR appears in Ea. (|59l) both in terms proportional to the 
graviton mass, and in the coefficient standing before H. The last fact is resulting 
from the proportionality of N to y^. Equations ([59| . so, do not reproduce the 
famous hypersurface deformation algebra [l8],[25, because functions N, do 
not coincide with parameters of this algebra. 

The substitution into Eqs. (|60)) expressions of the form (|55]) . corresponding 
to Poincare transformations, in the place of arbitrary functions a, a*, /3, /J-', 
results in Poincare algebra relations for the Dirac brackets: 

{P°,P^}d = 0, {P,,P,}d = 0, 
{P^M^'^lzj = 0, {P,,M^''}D^S,kPj~S,jPk, 

{pO,M*}z3 = -P\ {P„AP},5 = -<5„(F°-c"), 
{M^ M'^}d = 5kj{AP - c}) - 6m{M' - c^), {AP, AP}d = -AP^{61) 

Additive terms c° ~ tti'^/k J (fix andc* = m? / k J x^cfix in Pq and A'P, which 
do not depend on canonical variables and which are given by divergent spatial 
integrals play the role of central charges in Poincare algebra and correspond 
to the classical renormalization of the vacuum energy. They arise due to our 
intention to provide strictly zero energy for the empty Minkowskian spacetime. 
For this purpose a zero order over physical field term is introduced into the 
Hamiltonian (and also into the Lagrangian). As to linear over physical field 
terms in the Lagrangian and Hamiltonian, they appear only as total derivatives, 
and so do not contribute to the equations of motion. 

In principle, there is another opportunity: we can discard constant terms 
from the definitions of Poincare generators: 



P^ = -—I f^'d^X=- I n.d^X, 



K 
,2 



M 



ik 



^ J {x' f^'' - x'' f^') d^x = J {x'"Hi~x'nk)d^x, 



= -— I x'f^^d^x. (62) 



Then we get Poincare algebra without central charges: 

{P°,Pjz5 = 0, {P,,P,}d = 0, 

{P°,AP''}d = 0, {P,,A'P''}D^S,kPj~5,,Pk, 

{M'^,M'''}d = S,kM'' - SuAP'' + S.iAP'' - 5,kM'\ 

{P°,AP}z, = -P\ {P,,AP}d^-S^,P°, 

{M^M*^■},3 = 6k,AP-Sk^AP, {AP,AP}d^-AP'. (63) 
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But in this case the energy density for Minkowski space without matter wih 

2 

have a nonzero and positive value and the total energy will be infinite: 

= " V / ^^^'^''^ = V / ^^^^ 
We can conclude that it is natural to consider tensor 

2 2 / 

as a total energy-momentum tensor of gravitational field and matter. Let us 
add some reasons for this. The variational principle for RTG Lagrangian leads 
us to the following equations: 



where ^^atter energy-momentum tensor of matter 



1 TTi 

2^ 2 



2'^rnatten (66) 



2 SS 



^matter— 1 r • 1,0 V 

V~5 <^gyiv 

If we take a covariant derivative V^^, which is compatible with metric g^,^, of 
both sides of Eq. ([S5)) we obtain (compare with Ref. [H] Eq.(5.17)) 

^^total)-v ^matter: (68) 

where the semicolon denotes the covariant derivative compatible with the flat 
background metric. The matter equations of motion leads to equations 

V^T^'" 0, (69) 

so, if the gravitational field equations (|55|) are also valid, then we have 

(r^oD^.-O. (70) 

Then every Killing vector ^''(*), (i — 1,2,. . . ,10) of the Minkowski metric 

^I^i+ei:l. = o, (71) 

provides us with a conserved quantity 

\TZJi>dS,, (72) 

in correspondence with Eas. ((62| . The energy density (Eq.(57) or Eq.(62)) sign 
problem requires further study. Treated in the linear approximation the energy 
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density is an indefinite quadratic form, up to the spatial divergence. But only 
the full theory has an essential meaning. If some pathology in the form of 
negative energy flow generated by scalar component were present in the theory, 
then it should appear, for example, in the form of spherial waves radiation. 
But in publication [21] (see also Chapter 12 of Ref. [TS]) it was shown that 
the radiation of the gravitational field scalar component is abcent in spherically 
symmetric case and the external field is only static. Earlier analogous results 
were derived in Ref. [3S] 

6 Conclusion 

Let us summarize our results, formulate some expectations for the future and 
make some proposals. Up to now the massive gravity has not been considered in 
spacetime covariant Hamiltonian formalism (with arbitrary lapse and shift). As 
a result, Hamiltonians treated before did not provide freedom for the spacetime 
foliations. This work ensures us the Hamiltonian formalism flexibility which 
can be used to clear such fundamental issues as the causality in bimetric theory 
(see also bigravity [26]), as stability of the solutions, as positivity of energy 
and so on. In paricular, it is often required to compare predictions of massless 
and massive gravity for the same physical problem. As canonical variables and 
their brackets coincide, the main difference consists in four constraints which 
present in massless case and absent in massive one. Therefore, it would be rather 
interesting to compare the Cauhy problems with the same initial conditions for 
the two theories on rather long time interval, for example, compared with the 
Solar system age. We hope that the above mentioned problems are possible to 
study with modern numerical methods of computer calculations. 
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